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Abstract: 

The complete electroweak 0(a) corrections have been calculated for the charged- 
current four-fermion production processes e + e" — > u T r + ud/i - ^, and udsc. Here, 
technical details of this calculation are presented. These include the algebraic reduction 
of spinor chains to a few standard structures and the consistent implementation of the 
finite width of the W boson. To this end, a generalization of the complex-mass scheme to 
the one-loop level is proposed, and the practical application of this method is described. 
Finally, the effects of the complete 0{a) corrections to various differential cross sections 
of physical interest are discussed and compared to predictions based on the double-pole 
approximation, revealing that the latter approximation is not sufficient to fully exploit 
the potential of a future linear collider in an analysis of W-boson pairs at high energies. 
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1 Introduction 

W-pair production is an important process for testing the Electroweak Standard Model 
(SM). On the one hand, it allows to measure the mass of the W boson, a fundamental 
parameter of the SM, precisely. On the other hand, it is sensitive to the triple non-abelian 
gauge couplings and, in particular at high energies, allows to test the non-abelian structure 
of the SM accurately, owing to the delicate gauge cancellations in its lowest-order matrix 
elements. 

Experimentally, W-pair production has been studied intensively at LEP2 with quite 
high precision [[I]. The total cross section was measured from threshold up to a centre- 
of-mass (CM) energy of 207 GeV leading to a combined experimental accuracy of ~ 1%. 
While the W-boson mass Mw was determined from the threshold cross section with an 
error of ~ 200 MeV and by reconstructing the W bosons from their decay products within 
~ 40 MeV, deviations from the SM triple gauge-boson couplings were constrained within 
a few per cent. More precise experimental investigations of W-pair production will be 
possible at a future International e + e~ Linear Collider (ILC) [ El El E] . Owing to the 
high luminosity of such a collider, the accuracy of the cross section measurement will be 
at the per-mille level, and the precision of the W-mass determination is expected to be 
~ 10 MeV [Ej by direct reconstruction and ~ 7 MeV from a threshold scan of the total 
W-pair-production cross section [0 IS]- The higher energy of the ILC will enable much 
more precise measurements of the triple gauge-boson couplings, improving the sensitivity 
of LEP2 by more than an order of magnitude. 

Because of its theoretical importance, W-pair production found early interest. The 
lowest-order amplitudes for on-shell W-pair production were already considered at the end 
of the 1970's The electroweak corrections to on-shell W-pair production have been 
calculated by four different groups [ |H] shortly after and supplemented by hard photon 
radiation [ 9j. Already at that time, these calculations were at the forefront of the technical 
developments in higher-order calculations. Later, the structure of these corrections has 
been investigated by constructing improved Born approximations [ QUI and high-energy 
approximations [ ■ With the advent of LEP2 it became quickly clear that the decays 
of the W bosons into fermion pairs had to be included. The electroweak corrections to 
the on-shell W-boson decay were given in Ref. [H2|- Different types of programs (ranging 
from semianalytical codes to Monte Carlo generators) for lowest-order predictions for 
e + e~ — > 4/ were developed [[T3| (see also Refs. [ HH HE] and references therein) and 
subsequently supplemented by universal corrections thus reaching an accuracy of about 
2%. The universal corrections included running couplings, initial-state radiation, and 
also the effects of the Coulomb singularity for off-shell W-pair production [Ej. Since the 
accuracy of 2% was not sufficient for LEP2, the 0(a) corrections to e + e~ — > WW — > 4/ 
were calculated in the double-pole approximation (DPA), where only the leading terms 
in an expansion about the resonance poles of the two W-boson propagators were taken 
into account [ HH1 [THJ 1201 I2H 122] ■ These corrections were implemented into the event 
generators YFSWW [HH] and RacoonWW [ I2U 123 • A discussion of the remaining 
theoretical uncertainties of these calculations is presented in Refs. [El 121] • in view of 
the improved precision of the ILC, a further reduction of the uncertainties from missing 
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radiative corrections is necessary. This requires, in particular, the calculation of the full 
one-loop corrections for W-pair-mediated e + e~ — > 4/ processes. 

Such a calculation poses a number of theoretical challenges. 1 Neglecting diagrams in- 
volving couplings of Higgs bosons to light fermions, which are proportional to the fermion 
masses, already for the simplest final states e + e~ — ► v T r+ \i~v^ ud/x - ^, and udsc about 
1200 Feynman diagrams contribute (counting the contributions of the three fermion gen- 
erations in the loops only once), and for the most complicated final state u e e + e~u e , there 
are about 6000 diagrams. The large number and the complexity of the diagrams require 
to develop improved reduction algorithms, in order to keep the expressions manageable 
and to produce an efficient and numerically stable computer code. 

Because of the complicated multi-particle final state, the appearing loop integrals 
in general cannot be evaluated with standard methods. Using the Passarino-Veltman 
reduction to calculate the tensor integrals leads to serious numerical problems when the 
Gram determinants that appear in the denominators become small. This usually happens 
near the boundary of phase space but can also occur within phase space because of the 
indefinite Minkowski metric. Thus, in order to obtain numerically stable results, one has 
to devise and implement alternative methods for the calculation of the tensor integrals, 
at least, in the critical regions. 

The inclusion of the finite gauge-boson decay width constitutes a further important 
problem in the calculation of radiative corrections to W-pair-mediated four-fermion pro- 
duction. An appropriate description of resonances in perturbation theory requires a Dyson 
summation of self-energy insertions in the resonant propagators. It is well known that 
this procedure in general violates gauge invariance, i.e. destroys Slavnov-Taylor and Ward 
identities and disturbs the cancellation of gauge-parameter dependences, because differ- 
ent perturbative orders are mixed [ 121] • Several solutions have been described for 
lowest-order predictions. The early attempts have been summarized in Ref. [ 128] . and 
some of the schemes have been compared in Ref. [I2H]- More recent approaches include 
the "pole-scheme" [EDI EH, the "fermion-loop scheme" [ 123 HE E21 E3] , the use of effec- 
tive Lagrangians [EUEH] 2 , an d the "complex-mass scheme" (CMS) [123]- Apart from the 
pole expansions, none of these approaches has been elaborated beyond tree level so far. 
The pole scheme provides a gauge-invariant answer in terms of an expansion about the 
resonance, but is only applicable sufficiently far above the W-pair threshold. However, 
in the full calculation we are after a unified description that is valid everywhere in phase 
space, without any matching between different treatments for different regions. Some 
problems related to the finite width appearing in a calculation of radiative corrections to 
e + e~ — > yU~z^ud are illustrated in Ref. [EH]- Here we propose to solve these problems 
by using a generalization of the CMS, which was introduced in Ref. [|2B] for lowest-order 
calculations, to the one-loop level. 

In this paper we present details of our solutions to the above-mentioned problems. We 
discuss, in particular, our methods for the algebraic reduction of the one-loop amplitude 
and describe the use of the CMS at the one-loop level. In addition, we provide numerical 

^ome of the problems appearing in a first attempt of such a calculation were already described in 
Ref. [ESI- 

2 The recently proposed approach [ESI to describe unstable particles within an effective field theory is 
equivalent to a pole expansion. 
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results for the effects of the complete 0(a) corrections to the processes e + e~ — > u T r + jj,"^^, 
ud/z - ^, and udsc for various differential cross sections of physical interest. The effects 
of the complete corrections on the total cross section without cuts have already been 
discussed in Ref. [l3Tj. 

The paper is organized as follows: In Section El we fix our conventions and sketch the 
general strategy of our calculation. Our methods for the reduction of spinor chains to a 
few standard structures are described in Section EH Section |U is devoted to the description 
of the complex-mass scheme for loop calculations with complex masses. In Section we 
discuss the corrections for various distributions. We note that Sections EJ and El can be 
read independently. Section |H1 contains our conclusions. 

2 Strategy of the calculation 

The actual calculation builds upon the Racoon WW approach [ E] > where real- 
photonic corrections are based on full matrix elements and virtual corrections are treated 
in DPA. Real and virtual corrections are combined either using two-cutoff phase-space 
slicing or employing the dipole subtraction method as formulated in Ref. [ EH] for photon 
radiation. We also include leading-logarithmic initial-state radiation (ISR) beyond 0(a) 
in the structure-function approach (see Ref. [E] and references therein). The presented 
calculation only differs from RacoonWW in the treatment of the (IR and collinear 
finite part of the) virtual corrections. Therefore, we only describe the calculation of the 
complete 0(a) virtual corrections in the following. 

2.1 Notation and conventions 

We consider the process 

e + (p + ,cr + ) + e~(p_,<7_) -> f^h,^) + f 2 (k 2 ,a 2 ) + f 3 (k 3 ,a 3 ) + f 4 (k 4 ,a 4 ). (2.1) 

The arguments label the momenta p±, ki (i = 1,2,3,4) and helicities a±, <Ji = ±1/2 of 
the corresponding particles. We often use only the signs to denote the helicities. The 
particle momenta obey the mass-shell conditions p\ = m 2 and kf = mf. The masses of 
the external fermions are neglected whenever possible, i.e. everywhere but in the mass- 
singular logarithms. For later use, the following set of kinematical invariants is defined: 

s = {pi+p 2 ) 2 , ay = {ki + k,) 2 , t ±i = {p±-k t ) 2 , i, j = 1,2,3,4. (2.2) 

In this paper we consider only final states where f\ and f 3 are different fermions ex- 
cluding electrons and electron neutrinos; f 2 and f 4 are their isospin partners, respectively. 
This corresponds to the CC11 family in the classification of Ref. [E]- It represents the 
gauge-invariant subclass of general e + e~ — > 4/ processes that includes all diagrams with 
pairs of potentially resonant W bosons. In this class, the lowest-order and one-loop ampli- 
tudes vanish unless a = a~_ = — a + , <J\ = —<J 2 , and a 3 = — cr 4 . Moreover, the helicities of 
the outgoing fermions are fixed, = — 02,4 = — 1/2, owing to the left-handed coupling 
of the W bosons. In general this does not hold for other e + e~ — > 4/ final states and 
diagrams. We set the quark-mixing matrix to the unit matrix, but in the limit of small 
masses for the external fermions a non-trivial quark-mixing matrix can be easily taken 
into account by rescaling the cross sections for definite flavours accordingly. 
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Figure 1: Lowest-order diagrams for e + e — > /1/2/3/4 

The lowest-order cross section reads 

a = 1 |d$ 4 \( l + 2P-<j)(1 " 2P +°) l-^o^l 2 - (2-3) 
a=± h 

Here M.Q aius denotes the lowest-order matrix element, P± the polarization degrees of the 
beams, and d$4 the 4-particle phase-space volume element 

^=(np^»)(2^(p ++ P--|:%). (2.4) 

2.2 Survey of one-loop diagrams 

The virtual corrections receive contributions from self-energy, vertex, box, pentagon, 
and hexagon diagrams. In this section, we survey the diagrams contributing to the mass- 
less charged-current processes e + e~ — > /1/2/3/4, where fi and are different fermions, 
excluding electrons and electron neutrinos, and / 2 and / 4 their respective isospin partners. 

The contributions from self-energy and vertex corrections are obtained by inserting 
self-energies and vertex corrections in all propagators and vertices of the tree-level dia- 
grams shown in Figure^ When inserting a self-energy in a j/Z line one obtains 77 and 
ZZ self-energies as well as 7Z and Z7 mixing-energies. Since we neglect the masses of 
the external fermions, all diagrams that involve Higgs-boson couplings to these fermions 
obviously vanish. Nevertheless there remain diagrams containing contributions to the #7, 
HZ, and (f)W mixing energies and to the Hee and <ftff vertices with on-shell fermions, 
where denotes the would-be Goldstone boson corresponding to the W boson. One can 
easily verify that these contributions also vanish in the limit of vanishing masses for the 
external fermions. 

The diagrams for the appearing eej, eeZ, ev e W , ■yWW, and ZWW vertex functions 
are listed in Ref. [Ej, where the process e + e~ — > W + W~ was treated at one loop. The 
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diagrams for the other 7//, Zff, Wff vertex functions are simply obtained from the 
latter by obvious substitutions and adding some diagrams with internal Z bosons replaced 
by photons, which are absent for the ev e W vertex because of the vanishing charge of the 
neutrino. The diagrams for the gauge-boson and fermion self-energies can be found in 
Ref. [EH]- For all these vertex functions, the corresponding counterterm diagrams must 
be included. 

The generic contributions of the different vertex functions with more than three ex- 
ternal legs are shown in Figure 121 These are all ultraviolet (UV) finite. There are 40 
hexagon diagrams, 112 pentagon diagrams, and 227 (220) box diagrams in the conven- 
tional 't Hooft-Feynman gauge (background- field gauge [HO]). A set of hexagon diagrams 
is shown in Figure El and one set of pentagon diagrams for the Z/1/2/3/4 vertex function 
in Figure |U In both cases there are three further sets of diagrams that are obtained by 
reversing the fermion flow in one or both of the fermion lines of the outgoing fermions 
and by exchanging f\ <-> / 2 and/or / 3 <-> / 4 . Those for the 7/1/2/3/4 vertex function are 
simply obtained by replacing the external Z boson by a photon in the diagrams for the 
Z/1/2/3/4 vertex function. A set of diagrams for the eef\f 2 W vertex function is listed 
in Figure E] and a further set is obtained from those by reversing the fermion flow in the 
fermion chain of the outgoing fermions and by exchanging fi f 2 . The diagrams for 
the eefzf^W vertex function can be obtained from the latter by obvious substitutions. 
The box diagrams of the ZW f\f 2 vertex function are depicted in Figure El and those for 
the ZWfsfi vertex can again be obtained from those. The diagrams for the 7W/1/2 and 
7W/3/4 boxes are obtained by replacing the external Z boson by a photon and omitting 
the diagrams with internal H lines. The diagrams for the eeWW, /1/2/3/4, and ee/1/1 
boxes are compiled in Figures [3 El and 03 The ev e f\f 2 and ez/ e / 3 / 4 box functions are 
special cases of the j '1/2/3/4 box function. 

2.3 Calculational framework 

The amplitudes are generated with FeynArts, using the two independent versions 
1 and 3, as described in Refs. [ HI] and [E2]> respectively. The algebraic manipulations 
are performed using two independent in-house programs implemented in Mathematica, 
one of which builds upon FormCalc [US]. 

All contributions to the matrix elements involve three spinor chains, corresponding to 
the three different fermion-antifermion pairs, which are contracted with each other and 
with four-momenta in many different ways. There are 0(1O 3 ) different spinor structures 
to calculate for e + e~ — > /1/2/3/4 so that an algebraic reduction to a standard form which 
involves only very few standard chains is desirable. We have worked out algorithms 
that reduce all occurring spinor chains to O(10) standard structures, the standard matrix 
elements (SMEs), without introducing coefficients that lead to numerical problems. These 
algorithms are described in Section El 

It is convenient to separate the matrix elements into invariant coefficient functions 
F n , containing the loop integrals, and SMEs Ai n , containing all spinorial objects and the 
dependence on the helicities of the external particles [El]. After the reduction of the 
spinor structures, the amplitudes take the form 

M *^3 = ^Fr a3 ({s,s^t ±l })Mr a3 (P + ,P-,h,k 2 ,h,h). (2.5) 
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Figure 2: Contributions of vertex functions with at least four external legs to e + e 
/i 72/3/4 
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Figure 3: Ten hexagon diagrams for e + e~ — > /1/2/3/4. The remaining 30 hexagon dia- 
grams are obtained by reversing the fermion flow in one or both of the fermion lines of 
the outgoing fermions and by exchanging fi <-> / 2 and/or / 3 <-> / 4 . 




Figure 4: Six pentagon diagrams contributing to the Zfif 2 fzf& vertex function. The 
remaining 18 diagrams are obtained by reversing the fermion flow in one or both of the 
fermion lines of the outgoing fermions and by exchanging f\ <-> f 2 and/or / 3 <-> / 4 . 




7/^ ft 



l/Z h 



Figure 5: A set of 16 diagrams contributing to the eefif 2 W vertex function. The remain- 
ing 16 diagrams contributing to this vertex function are obtained from those by reversing 
the fermion flow in the fermion line of the outgoing fermions and by exchanging j\ <-> / 2 . 
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Figure 6: Diagrams contributing to the ZWA/2 vertex function 

The invariant functions F n are linear combinations of one-loop integrals with coefficients 
that depend on scalar kinematical variables, particle masses, and coupling factors. 

The calculation of the virtual corrections has been repeated using the background- 
field method [ED], where the individual contributions from self-energy, vertex, box, and 
pentagon corrections differ from their counterparts in the conventional formalism, apart 
from those that involve only fermion-gauge-boson couplings in the loops. The total one- 
loop corrections of the conventional and of the background-field approach were found to 
be in numerical agreement. 

The contribution of the virtual corrections to the cross section is given by 

5a virt = 1 J d$ 4 ]T ^(l + 2P_a)(l-2P + a)2Re{Mr-(M°-y}, (2.6) 

<T=±\ 

where A41 aia3 denotes the one-loop contributions to the helicity amplitudes. 

2.4 Calculation of loop integrals 

For the calculation of the loop integrals we use an improved version of the standard ap- 
proach. All coefficient integrals are algebraically reduced to a set of master integrals. The 
reduction formulas are implemented into Fortran codes, and the reduction is performed 



8 




Figure 7: Diagrams contributing to the eeWW vertex function 
fi l/ z h h W h h h h h 
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Figure 8: Diagrams contributing to the /1/2/3/4 vertex function 
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Figure 9: Diagrams contributing to the eefifi vertex function 
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numerically. Finally, the master integrals are either calculated from explicit formulas or 
from numerical integration as explained below. 

Let us first describe our standard approach that we use as long as it yields numeri- 
cally stable results. The 6-point integrals are directly expressed in terms of six 5-point 
functions, as described in Refs. 5-point integrals are written in terms of 

five 4-point functions following the method of Ref . [ HE] ■ The 3-point and 4-point tensor 
integrals are algebraically reduced to the (standard) scalar 1-point, 2-point, 3-point, and 
4-point functions, which are the master integrals, with the Passarino-Veltman algorithm 
[ I47j . Finally, the remaining scalar integrals are evaluated based on generalizations of the 
methods and results of Refs. [ EHJ IHH ED] • UV divergences are regulated dimensionally 
and IR divergences with an infinitesimal photon mass 3 m T 

The standard approach leads to serious numerical problems when the Gram deter- 
minants, which appear in the denominator of the Passarino-Veltman reduction, become 
small. This happens usually near the boundary of phase space but also occurs within 
phase space because of the indefinite Minkowski metric. Since we use Passarino-Veltman 
reduction only for 3-point and 4-point functions, the dangerous Gram determinants are 
those involving two or three momenta. Our reduction of 5- and 6-point functions, on the 
other hand, does not involve inverse Gram determinants composed of external momenta. 
Instead, it involves so-called (modified) Cayley determinants, the zeroes of which are re- 
lated to the Landau singularities of the (sub-)diagrams. We did not encounter numerical 
problems with these determinants. For the 2-point tensor coefficients numerically stable 
explicit results exist for arbitrary momenta [14*7]. 

For the regions where the Gram determinants become small, we have worked out two 
alternative calculational methods. The precise description of these will be given elsewhere. 
Here, we sketch only the basic strategy. One method makes use of expansions of the tensor 
coefficients about the limit of vanishing Gram determinants. Using this, again all tensor 
coefficients can be expressed in terms of the standard scalar 1-point, 2-point, 3-point, 
and 4-point functions. In the second, alternative method we evaluate a specific tensor 
coefficient, the integrand of which is logarithmic in Feynman parametrization, by numer- 
ical integration. Then the remaining coefficients as well as the standard scalar integral 
(numerator equal to one) are algebraically derived from this coefficient. This reduction 
again involves only inverse Cayley determinants, but no inverse Gram determinants. In 
this approach, the set of master integrals is not given by the standard scalar integrals any- 
more. For some specific 3-point integrals, where the Cayley determinant vanishes exactly, 
analytical results have been worked out that allow for a stable numerical evaluation. 

2.5 Checks on the calculation 

In order to prove the reliability of our results, we have performed a number of checks, 
as described in more detail in Ref. [EHj. We have checked the structure of the (UV, 
soft, and collinear) singularities, the matching between virtual and real corrections, and 
the gauge independence (by repeating the calculation in a different gauge). The most 
convincing check for ourselves is the fact that we worked out the whole calculation in 

3 At one loop the photon-mass regularization is equivalent to dimensional regularization (D = 4 — 
2e) with reference mass [i via the correspondence In(m^) «-> 1/e + ln(47r/i 2 ) — 7e as long as collinear 
singularities are regularized with fermion masses. 
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two independent ways, resulting in two independent computer codes the results of which 
are in good agreement. All algebraic manipulations, including the generation of Feynman 
diagrams, have been done using independent programs. In particular, one calculation uses 
the strategy described in Section I3~21 for the reduction of the spinor structures, the other 
the strategy presented in Section l3~3l For the calculation of the loop integrals we use 
the two independent in-house libraries which employ the different calculational methods 
sketched in Section 12.41 for the numerical stabilization. 



3 Algebraic reduction of spinor chains 

One-loop amplitudes for processes with six external fermions involve O(10 3 ) different 
spinor structures, which are products of three spinor chains. These structures are, how- 
ever, not independent, but can be related by using the Dirac algebra, the Dirac equation, 
momentum conservation, and relations that follow from the four- dimensionality of space- 
time. Of course, the latter relations can only be used after cancelling UV divergences, 
which are regularized dimensionally in our work. 

Below we describe two strategies that can be exploited to express any product of three 
spinor chains for a reaction involving six external massless fermions (and no external 
bosons) in terms of very few standard structures. Note that none of the following steps 
leads to Gram determinants in denominator factors, which potentially spoil the numerical 
stability 4 . For generic reasons we take all six fermion momenta pi as incoming (J2i=\Pi — 
0) and use the following shorthand notation for a spinor chain, 



.4 



± 



ab 



v a (p a ) Au±u b (p b ) , (3.1) 



where v a (p a ) and Ub{pb) are the usual spinors for the ( ant i-) fermions and 

w ± =^(l± 75 ) (3.2) 

the chirality projectors. Since we deal with massless fermions, the matrices A, which act 
in Dirac space, always consist of a product of an odd number of Dirac matrices 7^. The 
quantities that we want to simplify have the general structure 



v 1 (p 1 )Auj p u 2 (p2) x v 3 (p 3 )Bu a U4,(p 4: ) x v 5 (p 5 )Cu T u 6 (p 6 ) = ABC . (3.3) 

We start out by deriving and summarizing some basic relations that are used in the 
reduction. 

3.1 Basic relations 

3.1.1 Identities for Dirac matrices 

The four-dimensionality of space-time leads to the Chisholm identity 

-fyPjl = 0«0yr _ g^yP + gfhy* + ie a/3 7 5 T5T5) (3.4) 



'inverse Gram determinants, e.g., appear in the reduction proposed in Ref. [I51|. 
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which can be used to express products ^ a ^^"i of Dirac matrices in terms of single Dirac 
matrices dressed by 75 factors and totally antisymmetric tensors e a/3lS (we adopt the 
convention e 0123 = +1). 

If at least two Dirac chains are involved, this identity can be applied to shift factors 
of Dirac matrices from one chain to another without introducing e-tensors, provided the 
two chains are Lorentz contracted with each other: 

7Y7^ ® 1, ^ (/V - 9^1 a + <f V + ie^7,7 5 ) ® 7„ 

= ^ ® 7° - l a ® 7^ + <f V ® 7m - 7,7s ® (ie^^s) 7s 

^ 7^ ® 7° - 7" ® 7^ + <? Q/ V ® 7a. + 7°75 ® 7^75 - 7^75 ® 7°7s 

+ 3 Q/ V75 ® 7/,75 - 7 M 75 ® 7 / V7 / V (3.5) 

Here and in the following the symbol ® denotes the direct product in Dirac space, i.e. 
Dirac matrices are not sandwiched between Dirac spinors in this case. Relation (|3.5J) can 
be expressed in a very compact way after introducing the chirality projectors (|3.2|) and 
rearranging the order of some Dirac matrices, 

7 M 7 a 7 /3 u; ± <g> 7 M u; ± = ^l)± ® 7M7 /3 7 a ^±, 
7 a 7 /3 7 M cu ± ® 7 m cj ± = 7 M ^ ± <g> 7 /3 7 a 7A t ^±, 
7 M 7 a 7 /3 w ± <g> 7 m cl> t = 7 M u; ± <g> 7 a 7 /3 7 M ^ T , 

7 a 7 / Va; ± ® 7m^t = 7 M ^± ® 7/,7 a 7 / W ( 3 - 6 ) 

The second and the fourth relations are equivalent forms of the other two. 

Equations (J3.6|) lead to analogous relations for direct products of Dirac chains with 
more than one contraction. 5 For two contractions we get 

7^7 a 7^ ± ® J^J v uj ± = 4c/ a/ Vw± ® 7^±, 

7^7^ ® 7^7^ = <g> 7°w T . (3.7) 

The first of these relations is obtained by multiplying the second chain of the first relation 
of (J3.6|) by 7 7 7/3 from the right and renaming the indices /3 — > v and 7 — > (3 after some 
obvious algebra. The second relation is found analogously. 

Relations for direct products of Dirac chains with three contractions directly result 
from ()3.7j) upon contraction with g a p, 

YYY^± ® 7^7,7p^± = 167 M u; ± ® 7 M u;±, 

Yl u Y^± ® 7^7*7^ = 4 7 M ^± ® 7/^t- ( 3 - 8 ) 

Relations for more than three contractions could be derived from the above results recur- 
sively, but are not needed in our case. 



3 Such relations can also be found in Ref. [152]. 
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3.1.2 Identities for products of spinor chains 



Further interesting relations follow if we sandwich the Dirac matrices in the above 
relations between spinors and perform some simplifications. For instance, if we contract 
the second relation in ()3.7|) with pi } aP3,/3, then from the r.h.s. multiply the first chain of 
Dirac matrices by A and the second by B, where here A and B stand for the unit matrix 
1 or any product of an even number of Dirac matrices, and finally attach the external 
spinors, we obtain 



± r 



31 



34 



(3.9) 



On the l.h.s. we can eliminate j>\ and fa by the respective Dirac equations and obtain 
4:(pip 3 ) 7^ 'A ^ 7j,-B . Three similar relations can be derived from ()3.7|) with other 
suitable contractions. The resulting four relations read 



[Aj>; 

M 

Ah 



fcB 

iiB 
Bfc 



± 

34 
± 

34 

T 

34 
T 

34 



(p 2 p 3 ) [AY 

(pm) [fA 

(pm) \fA 
(p 2 p 4 ) \Af 



± r 



± 



JuB 

12 L ^ J 34 

± r i± 
£7 



B ltl 



34' 

=F 

34' 
=F 
34 ' 



(3.10) 



3.1.3 Decomposition of the metric tensor 

In four dimensions the metric tensor g^ u can be decomposed in terms of four indepen- 
dent orthonormal four- vectors nj, 



sT = E 9 ij < 



n ■ 



(3.11) 



i,j=0 



with orthonormal components rij^ = gij, where #y = = diag(l, — 1, — 1, — 1), as 
described in the appendices of Refs. [ESI Ell- The four four- vectors nf can be constructed 
from three linearly-independent momenta Pi,Pj,Pk- in the case of massless momenta, 
Pi — Pj — Pk = 0; they can be defined as 



n^ipi,Pj,Pk) = 1 1 (Pi+Pj 



ri2(Pi,Pj,Pk) 
n 3 {pi,Pj,Pk) 



ni(Pi,Pj,Pk) 



2 (PiPj) ipiPk) {PjPk) 
1 

2 (PiPj) {PiPk) {PjPk) 



2 (PiPj) 

{PjPk) Pi + (PiPk) Pj ~ (PiPj) Pk 



Pi ~P^ 



3 7 ' 



'Pi,aPj,f3Pk,j 



(3.12) 



The decomposition of the metric tensor (|3.11J) can be used to disconnect contractions of 
spinor chains or other objects. Note that the construction of the four vectors from the 
three independent vectors Pi, pj, pk, in particular the definition of 713, avoids an inverse 
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Gram determinant in (|3.11j) ; if we decomposed the metric into four totally independent 
momenta, the result would get the Gram determinant of those four momenta in the 
denominator. 



3.1-4 Identities involving e-tensors 



The fact that there is no totally antisymmetric tensor of rank 5 in four space-time 
dimensions leads to the Schouten identity 



e [a/3j5 9p]v 







(3.13) 



where [. . .] means antisymmetrization in a,[3,j,8, and \i. 

The product of totally antisymmetric tensors can be expressed as a determinant of 
metric tensors, 



gap 


g au 


gap 


g aa 




gf* 


g Pp 


9* 


9™ 


g^ 


glP 


9 la 


gSp 


g Su 


g Sp 


gSa 



(3.14) 



3.2 Strategy for reducing spinor structures 

In this section we describe one algorithm for the reduction of spinor chains. A second 
one is described in Section l3~3l The algorithm described here reduces all spinor chains 
to a minimal set. Although this algorithm is applicable more generally, we refer in the 
description to the massless charged- current 6f processes e + e~/i/ 2 / 3 /4 — > 0, where fi and 
/3 are different fermions excluding electrons and electron neutrinos and fa and f± their 
isospin partners. 



Step 1 Reduction of contractions between spinor chains 

As first step in the reduction, spinor chains contracted with other spinor chains or 
e-tensors are disconnected. 6 This is achieved by introducing a decomposition of the met- 
ric tensor (|3.11|) between contracted Dirac matrices or contractions between a totally 
antisymmetric tensor and a Dirac matrix: 

3 

1p ® 7 M = 9^ 1p ® lv = E 9 ij ii ® ii , 

i,j=0 

e^lp = ^ vp °9paT ™ ^ fi 9' J <>,,,».!■ (3-15) 

i,j=0 

The choice of momenta for a suitable decomposition of the metric tensor depends on 
the positions of the contracted matrices in the spinor chains. Preferably the momenta 
Pi, Pj,Pk entering (|3.12|) and thus (|3.11j) are selected in such a way that the Dirac equations 

6 The explicit e-tensors result from calculating the traces of Dirac matrices in diagrams with closed 
fermion loops. 
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v(Pi)fa = = $iu(pi) or the mass-shell condition fa\ = p\ = can be used most directly, i.e. 
without unnecessary anticommutations of Dirac matrices by means of the Dirac algebra, 
since additional terms come with it. Therefore, we count the number of anticommutations 
of the slashed momenta fa introduced by the decomposition of the metric tensor, with 
Dirac matrices of the chain that would be necessary in order to apply the Dirac equation 
or the mass-shell condition. We choose the three momenta with the fewest necessary 
ant icommut at ions . 7 

The following substeps are recursively applied until no contractions of spinor chains 
are left: 

• Disconnect a single Lorentz contraction between two spinor chains or between a 
totally antisymmetric tensor and a spinor chain by (|3.15|) . 

• Use the Dirac algebra together with the Dirac equation or the mass-shell condition 
fa = in order to shorten spinor chains. 



Replace contractions between totally antisymmetric tensors and n^: 

QUI l 



[a /3 7] 
Pi PjPk, 



'2 (pipj) (piPk) (pjPk) 
where [. . .] means antisymmetrization in a, (3, and 7. 
• Eliminate ^3 in spinor chains by using 

jpj i [Afafa - (piPj) fa + (piPk) fa - (pjPk) fa] 75 



(3.16) 



MPhPjiPk) 



2 {PiPj) {PiPk) (PjPk) 



(3.17) 



After this step, the (disconnected) spinor chains we are left with are of the form \j> 



3 lab 



with j 7^ a,b and [^ifafa] ab with pairwise different i,j,k,a,b for processes involving six 
external fermions. Since there are many different types of multiple contractions of spinor 
chains, we can only illustrate this reduction step, and particularly the first substep, in 
three representative examples: 



ab 



7m 

rem 



cd 
1 

2 



3 

i,j=0 



1 r 



(PaPc 



fal5 



ab 



fa(Pa,Pb,Pc 

fa 
1 



ab 



cd + (p b p c 



ij(Pa,Pb,Pc) 



cd 



cd 



(PaPc 



fal5 



cd 



(PbPc) 



fal5 



cd 



More precisely: If the metric tensor between two spinor chains is eliminated, we count the necessary 
anticommutations for each momentum in both spinor chains, say Ni and N2 with N\ < N2, and set -/V2 
to infinity if the momentum does not appear in one of the chains. Then we choose the momenta with 
smallest N\ — I/N2. If, on the other hand, the metric tensor between a spinor chain and an e-tensor is 
replaced, we proceed analogously, but set the number of necessary anticommutations to zero for those 
momenta that are already contracted with the e-tensor and to infinity otherwise. 
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r 



12 



fa 



fa 7^ lv 



ab \(PaPc 



fa . $ct,-t + 7 

J cd [PbPc 



fa J<r,r 
cd 



(3.18) 



7 



56 



9 13 \ii(p2,pi,P3)\ l2 \faij(p^PuPz)i> 



i,j=0 



3d 



7 



56 



(PiPs) 

i3^T5l 1 

(PiPs) 
isTTl 1 



12 



fa 1n j*l ^2 7^ 



31 



7 



56 



7„ 



^] 12 £ ^ [^1^2^i(p2,Pd,P5)] 34 [ij{P2,P4,P5 
i,j=0 

fa]^ \fafafa] OA \fa 



56 



(3.19) 



12 



ie^pi^^p^ 9 %3 ni,i*(pi>P2,P3) [^j(Pi,P2,Ps) 

i,j=0 



12 



iSTSl 1 
FTT71 2(p 2 p 3 ) 

1 



fa\ 12 ( (pipa) (psPd) - (pm){p2Pi) + (pm)(P2P3) 



-i e a/37<5 pi,aP 2i /3P3 i7 Pd,«5 



2 (P2P3) 



^3 



12 ^123d J 



(3.20) 



where we introduced the abbreviation 



Kj ki = a (PiPj) (PkPi) + b (piPk) (pjPi) + c (pipi) (pjp k ) -die ijk i (3.21) 

with eijki = e a j3ryspfpjplpi and upper index combinations ( — h +±), (H h±), and 

(+ + -±) for (abed). 

Step 2 Reduction of a spinor chain to standard form 

Since all spinor structures are disconnected after the first step, we can focus on single 
spinor chains in this step, i.e. spinor chains of the form [B]° b where all Dirac matrices 
in B are contracted with momenta. These spinor chains can be reduced to the standard 
form \f^ b with a freely chosen momentum p = p n , p n ^ p a ,Pb, by recursively applying the 
following substeps: 

• For p m ^ p a ,Pb,P, replace j> m by 



fan = Pm,n 9 l ° n% ij 

i,3=0 



(3.22) 



Eliminate ifi^ via ()3.17|) . 



Use the Dirac algebra together with the Dirac equation or the mass-shell condition 
fa = in order to shorten the spinor chain. 
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Since for our case only two different types of spinor chains are left after the first step, we 
can demonstrate the reduction procedure of the two cases in full detail. The first case is 
reduced as: 



in 



± l3~2"2l 

ab 



i,j=0 

(233 {PaPn) {PbPm) ~ {PaPb) {PnPm) + {PaPm) {PnPb) ± 1 tanbm 



1 



2 (PaPn) (PbPr. 



a nbm 



2 (PaPn) (PbPn) 
fa] 



fa 



± 

ab 



ab 



The second case can be reduced to the first follows: 



(3.23) 



i>m i>\ fa 



-j- ^^^^^ ^ -j -|- 

Pl,H 9* 3 U i(Pm,Pk,Pb) imij(Pm,Pk,Pb)ik 



ab 



i,j=0 



ab 



EI3 (PbPm) (PlPk) - (PbPl) (PmPk) + (PbPk) (PmPl) 



— 1 



tbmlk 



2 (PbP m )(PbPk) 



i> m fafa 



2 (PbPm)(PbPk) 
(PbPm) {PlPk) ~ {PbPl) (PmPk) + (PbPk) (PmPl) =F '^^bmlk 



(PbPr. 



fa, 



l3~2Tl 



{PbPr. 



4+-+± 

^bmlk 



fa, 



(3.24) 



Thus, we end up with spinor chains in standard form and prefactors S containing only 
scalar products of external momenta (piPj), Cijkii an d the abbreviations (j3.21)) . Note that 
()3.24|) could be applied recursively to longer spinor chains. 

Thus, for processes with six external fermions, all Dirac structures occurring in one- 
loop amplitudes, i.e. all structures of the form ()3.3|) . can be brought to the form 



S 



12 



fa , fa ,„ i> 



31 



5(5 



(3.25) 



where we have chosen the momenta p 3 , pi, and p\ in the spinor chains by convention. 
Accordingly, the SMEs can be chosen as 



M 



par 



fa fa fa 

ro i 12 L rl J34 l r 



56 



(3.26) 



and there is only one SME for each helicity combination of the external fermions. For 
purely W-mediated charged- current processes, there are only two non-vanishing helicity 
combinations and thus only two different SMEs. 

Step 3 Simplification of scalar factors 

After the reduction steps above, each spinor structure has the form ()3.26|) with prefactors 
S containing products of scalars products (piPj) in the denominator, and possibly poly- 
nomials of scalar products, of e^^, and of A^lf in the numerator. These prefactors can 
be simplified further by means of relations we describe now. 
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The quantity A^ k f denned in (j3.21|) transforms under exchange of two momenta, 
corresponding to two of the indices (ijkl), as 

Aabc± Aacb^p Aacb^p Acba^p Acba^p Abac^p Abac^p jo nn\ 

A ijkl ~ A jikl ~ A ijlk — ^kjil ~ ^ilkj —^Ijki —^ikjl > 

and is therefore invariant under exchange of two distinct pairs of momenta, 

Aabc± Aabc± Aabc± <\abc± /o oo\ 

A ijkl — A jilk ~ A klij —^Ikji • {O.ZO) 

Owing to these relations, any A^^,, where (i'j'k'V) is an arbitrary permutation of (ij k I), 
can be transformed into one of the six elements 

Aijti ) A'ljM , or A~lj kl . (3.29) 

In the following, we use only these independent quantities (J3.29)) . i.e. for each set of indices 
(ij kl) of A1 b j C k f we define a standard order. 

The identity for products of totally antisymmetric tensors (|3.14j) leads to relations 
among the A^ k f. Relations for products of Aft-fif with the same momenta are 

A ~ijkl ' A ~ijM = 4 (PiPk)(PjPl)(PiPl)(PjPk), 

A + ijkt A+ ijkf = 4(piPj)(pkPi) (j>iPi)(pjPk) , 

A+ ijkt A YfkJ = 4 (piPj)(pkPl) (PiPk)(PjPi) , (3.30) 



and 



AtiktMtkt = -2{p lPj ){pm) A ^ + M 

+ + ± A++-± _ Of„„.\f„„.\ A+- + T 



A ~iUr A+ ijkf = -z{PrPk)(pm) A T 3 + k z i 

A T 3 tt Atjtt = - 2 MM A+ i+k? ■ ( 3 - 31 ) 
Note that in these and the following formulas double Latin indices are not summed. In 
()3.30|) and (|3.31|) the second and third lines are obtained from the first by the substitutions 
(j ^ k) and (i <-> k), respectively, and subsequent transformation to the six elements 
(I3~2H . 

Relations for products of A^^f 1 which differ in one momentum read 

a-++± — i A~ ++± A~ ++T = 1 A ++ ~ ± A ++ ~ T f3 32") 

^ijkl 2( Pi p m )(p jP m) ijkm^ijlm 2(p k p m )(pip m ) ^iklm ^jklrm yu.u^j 

and 

A+-+±A++-± _ (PiPj) A~++± = (P'Pi) _ (PiPi) A-++T A~++± 

^ijkl^ijkm ' {pip m ) i 3 km^ 1 i klm (pip m ) « k 1 m 3 klm (PjPl) ' 3 kl ^ J klm ■ 

(3.33) 

Two further relations can be obtained from each of the equations in 1)3.32)1 and (j3.33j) 
via the substitutions (j <-> k) and (i <-> k). Three further sets of relations can finally be 
constructed by substituting (I <-> m) in all these relations derived from ()3.33|) . 
Step 3 consists of two parts: 
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and linear relations like A^-^f 1 



First, we try to eliminate sums containing (p-iPj) (PkPi), tijkh an d products 

thereof. To this end, we use the Schouten identity (|3.13|) . the relations (|3.3U|) - (j3.33|) . 

2 (PiPj) (PkPi) ~ A + /fcT or (PiPj) (PkPi) ±ie ijkt = 
(PiPk) (PjPi) + (PiPi) (PjPk) that follow from the definition of the A^^i 
()3.21j) . All these relations are applied recursively to parts of S in ()3.2(ij) as long as 
they lead to simplifications. 

When starting from structures of the form ()3.3|) . we succeeded in this way to elimi- 
nate all sums, and the resulting S involve only products of A^^f and scalar products 
in the numerator, and products of scalar products in the denominator. 

• Second, the products of A^lf are simplified further and brought into a standard 
form by using the relations ()3.30|) - (j3.33l) . 

Typical results are 



12 



7m 



12 



12 





34 . 


7^7 P 7 K 


34 . 



12 

7^7p7« 





12 




34 


>. 




>. 


+ 

12 




34 


>. 



34 

56 
56 

56 
56 



(PlPs) 

0, 
0. 



2A\ 



234 



34 V (PlP4) {P2PZ, 



^ /1 1345 /1 1346 



(P1P3) (P1P4) (piPs) (PlPe) (P3P4) 
a/1 1234 /1 1356 

(pips) (pm) (pm) (pm) (pm) 

/1 1235 



56 2 (pjp 3 ) (pip 5 ) (p 2 p 3 ) 



12 



r i 34 1/ 1 



56 



^ _ j^l Q . 2* 



12 



56 



31 



56' 

(3.34) 



The entire reduction algorithm, described in the three steps above, reduces the spinor 
chains of the considered processes to the SMEs (J3~27)|) and about 35 different Af^f . When 
inserting the results into the amplitudes, further simplifications between contributions of 
different spinor structures can be performed owing to the simple structure of the SMEs. 
This speeds up the calculations considerably. 

The spinor structures could alternatively be evaluated directly in the Weyl-van der 
Waerden spinor formalism, as described in some detail in the appendix. Although most of 
the relations given in this section can be easily derived using the spinor formalism, many 
simplifications that are based on four-momenta are harder to perform in that approach. 



3.3 Alternative strategy for reducing Dirac structures 

In this section we describe an alternative strategy for the reduction of Dirac struc- 
tures. In this approach we take care that no quantities are introduced in the denominator 
that can lead to numerical problems, like scalar products or Gram determinants. As a 
consequence, the Dirac structures are not reduced to the minimal set, but to a small set 
of standard structures. Moreover, the coefficients in front of these standard structures 



19 



only involve Lorentz products, but no contractions with e-tensors, facilitating further 
simplifications in amplitudes. 8 

Step 1 Reduction of multiple contractions between spinor chains 



Multiple Lorentz contractions between two spinor chains can be reduced to single con- 
tractions using ()3.7|) and ()3.8|) . A typical example is given by 



AYiiY i2 Birffy. 



± K7\ 

34 



A(p iPj ) AY 12 B % 



± r 



34 



(3.35) 



This procedure leaves only products of spinor chains with (i) no contractions, (ii) exactly 
one contraction between two spinor chains, or (iii) the case in which one spinor chain is 
singly contracted with each of the two other chains. 

Case (iii) can be reduced to (i) and (ii) as follows. Each spinor chain of type (iii) can 
be brought to the form [AYi>iY]\2 [-^7^34 [^TfJIe or ^° one °f ^he ana l°g° us forms with 
external fermions interchanged. The transformation to this standard form produces only 
terms of types (i) and (ii). The factor fa in the first chain is always present, because the 
number of Dirac matrices in the chain is at least three. We can assume that Pi is the 
momentum of a spinor belonging to the other chains, i.e. i = 3,4,5, or 6. Otherwise we 
could eliminate it with the Dirac equation after anticommuting it to the left or right in 
the spinor chain [. . .] 19 . The trick to get rid of the two Lorentz contractions is to apply 
(13. 6|) in order to shift the factor ^ to the spinor chain where it can be eliminated with 
its Dirac equation. This means, if i = 3 or 4, fa is shifted into [. . .] 34 ; if i = 5 or 6, it is 



shifted into 



-• • -J56 



. We give a simple example: 



12 

AY 
AY 



31 



56 



AY\ 12 [BfaYlv 



31 



c 7 „ 



56 



12 
± r 



B(2p^ - 2p 4l/ i7 

T 



12 



Bl, 34 Cfa 



56 



- 2 



T 
34 

Ah 



± 

12 



56 



by c % 



31 



56 



(3.36) 



Step 2 Reduction of single contractions to standard form 



After the first step, the only contractions between spinor chains are of the form 
[-47^] 12 [-^7^34 [^56 ( or w ith fermions interchanged), where the matrices A, B, C are 
products of slashed momenta. The matrices A and B can be systematically eliminated 
as follows. 

We reduce the number of slashed momenta in A until A — 1. If A = 1, there is 
nothing to do. Otherwise A is of the form A'tfify with A' containing two slashed momenta 
less than A. If either % or j is 3 or 4, we shift j>^j>j to the spinor chain [. . .] 34 with 
the help of (J3.6)) and eliminate ^ 3 or ^ 4 with its Dirac equation. This leaves the cases 
i>ii>j = ^5^6 or ^6^5 (since p'ip'i = 0), which can be reduced by momentum conservation, 
e.g., = ^5 (— j>\ — j>i — j>z — i>i). The terms ^3 and j>\ are again shifted to [. . 
while j>\ and ^2 can be eliminated with their Dirac equations in [. . \xi directly. The whole 

8 In this approach all explicit e-tensors are eliminated with the Chisholm identity Ij3.4|l at the beginning 
of the reduction. 
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procedure reduces the number of slashed momenta in [Ay^]^ by two and can be repeated 



until A = 1. 

The same procedure can be applied to [-87^34 until B = 1. Finally, all products 
of spinor chains containing a Lorentz contraction between two chains are of the form 
[7 M ]i2 [7/J.34 [^56 ( or with fermions interchanged). Again, we give an illustrative example: 



fafaY 



fafaln 
..1 ± 
12 

2(P4Ps) 



C 



56 



r 



fafai^fafa 



31 



r 



± 



7 



fafain 



c 

34 



56 

c 



56 



= -2(p 4 P5) 
^ -2(p 4 Ps) { 
= -4(p 4 p 5 ) { 

+ (P4Pe) 



± 
12 

± 

12 



7' 



31 
A' 



7p 



± 
12 



31 



+ 



7" 



3d 
± 

12 



c 



56 



(fa + fa) fall] J 4 } 



c 



56 



± 






12 




34 



± r 



12 



7u 



T 
34 



^4] \fa] T } 
r J 12 L ru J 34 J 



C 



56 



(3.37) 



Step 3 Reduction of products of slashed momenta 



First we write all products of slashed momenta in a standard form. After the preced- 
ing steps, we are left with spinor chains of the form [A] p 12 [7^34 [7/J5 6 , [7 M ]i2 [-^34 [in] mi 
[7^12 [7J34 [^56' or [A]u [B]m [C]l, 6 - None of the A, B, C contains open Dirac indices, 
i.e. they are products of slashed momenta and chirality projectors. With the help of 
Dirac algebra and of momentum conservation the spinor chains can be brought into the 
form of a linear combination of terms \j)if ab and [fafafa]^- F° r [^]i2> e -&-; we substitute 
p e = —pi — p 2 — P3 — Pi — P5, eliminate all fa and fa with their Dirac equations by an- 
ticommuting them to the left or right and reorder the remaining fa terms according to 
increasing i. Owing to the relation fa = p 2 , this leaves only [fa p 12 with i = 3,4,5 and 
[fafafa]i 2 - In the same way we reduce the other uncontracted spinor chains to [fa] p Z4 , 

[fafafaYw DMb6> and [fafafa]l& with j = 5 ' 6 ' 1 and k = 2 > 3 - 

We now eliminate all chains of the form [fafafa] p ab [fa^>mfa]cd making use of the 
Chisholm identity ()3.4|) which can be written as 

fafafa = E ijk + G ijk (3.38) 

with the shorthands 



E 



ijk 



k al3l5 pi, a Pj,f3Pk r/ l8l5, 

(PiPj)fa - {PiPk)fa + (PjPk)fa- 



(3.39) 



Thus, we can write 



fa.fafa , fai> m fa 



a 
cd 



E. 



ijk 



p r 



+ 



Gin 



Elmn 
P 



+ 



ab 



cd 

i>ii> m fa 



fafap 1 , 



cd 



ab 
Gijk 



Glmn 
P 



ab 



G 



a 
cd 

Iran 



cd 



(3.40) 
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In the first term on the r.h.s. we can express the product of e-tensors, 



E„ 



E, 



Imn 



cd 



Pi,aPj,l3Pk,"/Pl,a>Pm,P>Pn,-y> 



7575 



75' 75 



cd 



(3.41) 



in terms of ordinary Lorentz products because of 1)3. 14jl . This means, after ex- 
ploiting (J3.40)) products of slashed momenta only occur in terms of the forms: (i) 

\fi3fafaY12 M34 Wse, (ii) \fizfafa\i2 [faYzA ^56' and sinmar terms with products of slashed 
momenta in the other spinor chains. 

Terms of type (ii) can always be simplified with (j3.1U|) . If one of the indices j, k is 1 
or 2, then fap 1 ^fa can always be reordered such that one of the relations in (J3.1U)) applies 
(if necessary momentum conservation has to be used). We illustrate this by a simple 
example: 



fahfe ln fa , A i>. 



12 
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34 



56 
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(3.42) 



If neither j or fc has the value 1 or 2 in terms of type (ii), then one of the relations in 
()3.10|) applies to the product [$7] 34 [fa]l e (if necessary momentum conservation has to be 
used). In this case terms of type (ii) transform into terms of type (i) and terms without 
products of slashed momenta. 

Finally, we have to reduce terms of type (i), which turns out to be a tedious task. We 
illustrate the basic steps for the following example, 
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fai>a v iAfa +fa+fa+fa) 
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31 



7 t 
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(3.43) 



At first sight, we have made the expression more complicated by using the Dirac equation 
and some rearrangements for fa, which was then eliminated by momentum conservation 
in the second equality. The last expression in 1)3.43)1 can be simplified by anticommuting 
the terms in (fa + fa + fa + fa) in such a way that the Dirac equations for fa and fa 
apply and that the terms fa and fa drop out via fa fa = 0. The various terms that survive 
this procedure can all be further reduced as described in the previous steps. After some 
tedious algebra, we get 
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Note, however, that the explicit form of this result is not unique, i.e. it might look different 
if other algebraic manipulations had been made. 



Step 4 Reducing products of the form [fa] p 2 [fa] 34 [A] 5, 



56 



Obviously, momentum conservation can be used such that the indices i,j, k take only 
the following values, i = 3,4,5, j = 5,6,1, and k = 1,2,3. Equations ()3.10|) provide 
2x3x3 = 18 relations per chirality which can be used to reduce the number of these 27 
(ijk) values to 9. 

There are, however, still some relations among the structures [A]^ [^j] 34 IfaJse and 
[fa]i2 [7I34 [7^56 j e ^ c -' wn i cn can be exploited to eliminate some (ijk) values. Such rela- 
tions can be deduced as in the following example, 
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This expresses the last term in terms of [A] 12 [7 M ] 34 [7m] 56 > e ^ c - ^ n total we find 7 such 
relations per chirality, so that the number of structures [A] 12 [fa] 34 [A] 56 can he reduced 
to 2 combinations (ijk) per chirality configuration (par). 

Step 5 Reducing products of the form [fa] p 2 [7^] 34 bvJse; e ^ c - 

As in the previous step, we use momentum conservation to constrain the indices i,j, k 
in the structures \^] p 2 [Yf u [7^, [7% [7^, and [7% to the values 

z = 3,4,5, j = 5,6,1, and k = 1,2,3. One out of these 9 resulting structures (per 
chirality) is, however, redundant and can be easily eliminated as follows, 
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Relations for the other chirality configurations can be obtained in the same way. 

We note that there are further relations among the spinor structures considered in 
this step. The remaining relations, however, involve coefficients with Lorentz products, 
so that their use to further eliminate some structures would lead to Lorentz products in 
the denominator, which we want to avoid. Nevertheless we give two such relations for 
illustration: 
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Summary: Final set of spinor structures 

After applying the steps described above, all Dirac structures occurring in one-loop am- 
plitudes for processes with six external fermions can be expressed as linear combinations 
of the following SMEs with coefficients that do not involve any denominators, 
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4 The complex mass scheme at one loop 



The description of resonances in perturbation theory requires at least a partial Dyson 
summation of self-energy insertions. This leads to a mixing of perturbative orders and, 
if done carelessly, can easily jeopardise gauge invariance [|2Ell2Zj- Therefore, the proper 
introduction of finite-width effects is a non-trivial problem. While several solutions have 
been described for lowest-order predictions [1231 123 12HI 12^1 EOI EH E21 E31 EH , no viable, 
universally valid scheme exists so far for a consistent evaluation of radiative corrections in 
the presence of resonances. A pole expansion [ |2Hl 1301 EH ES] provides a gauge-invariant 
answer and is applicable to radiative corrections, but restricts the validity of the result to 
the resonance region only and is not reliable in threshold regions. In our calculation we 
want to cover both the threshold region, where the pole approximation is not applicable, 
and the continuum above threshold, where threshold expansions are not valid. Moreover, 
the calculation should be valid both for resonant and non-resonant regions in phase space. 
In other words, we are after a unified description that is applicable in the complete phase 
space and does not require any matching between different treatments for different regions. 

Such a description is provided by the "complex- mass scheme" (CMS), which was 
introduced in Ref. [122] for lowest-order calculations. In this approach the W- and Z- 
boson masses are consistently considered as complex quantities, defined as the locations 
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of the poles in the complex k 2 plane of the corresponding propagators with momentum 
k. Gauge invariance is preserved if the complex masses are introduced everywhere in the 
Feynman rules, in particular in the definition of the weak mixing angle, 

cos 2 #w = 4 = l-4 = 4> (4- 1 ) 

which is derived from the ratio of the complex mass squares of the gauge bosons, 9 

/4 = Ml - iM w r w , /4 = Ml - iM z r z . (4.2) 

The (algebraic) relations, such as Ward identities, that follow from gauge invariance re- 
main valid, because the gauge-boson masses are modified only by an analytic continuation. 
As a consequence unitarity cancellations are respected, and the amplitudes have a decent 
high-energy behaviour. 

While necessary in the resonant propagators, the consistent introduction of complex 
gauge-boson masses introduces spurious terms in other places, as e.g. in the weak mixing 
angle (|4.1j) . When using the CMS at tree level, which amounts to replacing the real 
gauge-boson masses by the complex masses (J4.2)) and the weak mixing angle by (J4.1|) in 
tree-level amplitudes, the spurious terms are of order 0(Tw/M-w) = 0(a) relative to the 
lowest-order term (both in resonant and non- resonant regions). 

Here we propose a generalization of the CMS to higher orders. The complex masses 
are introduced directly at the level of the Lagrangian by splitting the bare masses into 
complex renormalized masses and complex counterterms. This scheme has the following 
properties: 

• From the Lagrangian we obtain Feynman rules with complex masses and countert- 
erms with which we can perform perturbative calculations as usual. Since we do not 
change the theory at all, but only rearrange its perturbative expansion, no double 
counting of terms occurs. 

• For each unstable particle mass, we add and subtract the same imaginary part 
in the Lagrangian. One of these terms provides the imaginary part for the mass 
parameter and becomes part of the free propagator, while the other becomes part 
of a counterterm vertex. The first term is, thus, resummed but the second is not. 
Independently of the imaginary part that is added and subtracted, this procedure 
does not spoil the algebraic relations that govern gauge invariance, and unitarity 
cancellations are exactly respected. In practice, this means that we can insert values 
for the gauge-boson widths that are not directly related to the one-loop order to 
which the corrections for the process are calculated. We could even go beyond one 
loop in the calculation of the widths or take an empirical value. 

9 While it is generally accepted that the mass and width of unstable particles are related to the pole 
of the propagator in the complex plane, this does not define the mass and the width separately. This 
arbitrariness is discussed in detail in Ref. [I55|. where also a definition of the mass and width is proposed 
such that the width is given by the inverse lifetime. We have chosen the popular definition Q4.2p . but the 
complex renormalization scheme is applicable to other definitions as well. 
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• Performing an 0{a) calculation in the CMS yields 0(a) accuracy everywhere in 
phase space provided the width that enters in the resonant propagators via the 
complex mass is calculated including at least 0(a) corrections. This is evident away 
from the resonances, where one could expand in terms of the width, thus recovering 
the usual perturbative expansion. In the resonance region, where the resonant 
contributions dominate, both the prefactors of the resonant propagators and the 
resonant propagators themselves are taken into account in 0(a) and our results 
differ by 0(a 2 ) terms from a leading pole approximation where this is applicable. 
Thus, any spurious terms are of order 0(a 2 ). 

Introducing complex masses and couplings seems to violate unitarity. Obviously, the 
Cutkosky cutting equations [ EH] are no longer valid, and unitarity cannot simply be 
proven order by order anymore. However, since we do not modify the bare Lagrangian, 
the unitarity-violating terms are of higher order, i.e. of 0(a 2 ) in an 0(a) calculation. 
Moreover, this unitarity violation cannot be enhanced, because all Ward identities are 
exactly preserved. In this respect one should also mention that unstable particles should 
be excluded as external states and only the S"-matrix connecting stable particle states 
needs to be unitary, as has already been pointed out by Veltman in the sixties [ EZ] • 
Of course, before the described CMS can be viewed as a rigorous procedure to define a 
renormalized quantum field theory it has to be clarified whether one can directly prove 
unitarity order by order in this formalism. In particular, it is an interesting question 
whether one can construct modified cutting equations in the CMS. 

(i) Complex renormalization - 't Hooft-Feynman gauge 

The consistent introduction of complex masses in loop calculations necessitates the 
formulation of an appropriate renormalization prescription. To this end, we generalize 
the on-shell renormalization scheme formulated in Refs. [ HHl EH EH at the one-loop level 
in a straight-forward way. A generalization to higher orders should be possible, but this 
is beyond the scope of this paper. 

Following the conventions of Ref. [Elj, the renormalized transverse (T) gauge-boson 
self-energies read 

£¥(k 2 ) = Y%(k 2 ) - 5/4 + (k 2 - n 2 w )5Z w , 
£ zz (k 2 ) = Z zz (k 2 ) - 5fi 2 + (k 2 - ^ 2 )5Z ZZ , 
t AA (k 2 ) = E AA {k 2 ) + k 2 5Z AA , 

t AZ {k 2 ) = H AZ {k 2 ) + k 2 h)Z AZ + (k 2 - A) l -5Z ZA) (4.3) 

where A denotes the photon field, and the hat indicates renormalized self-energies. Com- 
pared to Ref. [ 03] , the renormalized on-shell masses and mass counterterms are replaced 
by the renormalized complex masses /iw and \iz everywhere, i.e. also within the self- 
energies. We denote the field renormalization constants in the CMS by calligraphic let- 
ters. The complex renormalized masses and mass counterterms result from a splitting of 
the real bare masses squared, 10 

^w,o = /4 + <^w> M z 2 o = t4 + <5/4 (4-4) 

10 Similar ideas were proposed in Ref. 
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where here and in the following bare quantities are indicated by a subscript 0. Similarly, 
splitting the bare fields in complex field renormalization constants and renormalized fields 



Wt = (1 + 16Z W )W±, 

1 + \bZ zz \8Z ZA 
\5Zaz 1 + \5Zaa j \ 




A 



(4.5) 



implies that the bare and renormalized fields have different phases. Thus, for instance, the 
renormalized Z-boson field becomes complex, while the corresponding bare field is real. As 
a consequence, the renormalized Lagrangian, i.e. the Lagrangian in terms of renormalized 
fields without counterterms, is not hermitian, but the total Lagrangian (which is equal to 
the bare Lagrangian) of course is. 

In order to fix the counterterms, we generalize the renormalization conditions of the 
complete on-shell scheme [ EU EH| and require 





= o, 


X#Vz) 


= o, 






= 0, 




= o, 






= 0, 




= 0, 


e^(o) = o 



(4.6) 

where the prime denotes differentiation with respect to the argument. The conditions 
(|4.6|) . in particular the first two, fix the mass counterterms in such a way that the renor- 
malized mass is equal to the location of the propagator pole in the complex plane. This 
is a gauge- invariant quantity, as pointed out and shown in Refs. [EDJEDl- The last five 
renormalization conditions in (|4.6|) fix the field renormalization constants. Note that 
the field renormalization constants of the gauge-boson fields exactly drop out in all S- 
matrix elements that do not involve external gauge bosons, but allow to render all vertex 
functions finite. This generally holds for all field renormalization constants of unstable 
particles as long as one does not consider S-matrix elements for external unstable parti- 
cles. Unlike in Refs. [ EU EH] , we did not take real parts in the renormalization conditions 
([4.6)1 . and thus not only the mass renormalization constants but also the field renormal- 
ization constants become in general complex. This ansatz is supported by the fact that 
the imaginary part of one-loop scattering amplitudes involving unstable external particles 
becomes gauge dependent if the imaginary parts of the counterterms are not included [ 
l6*T| . For the definition of the renormalized mass and width this scheme is exactly the 
one described in Appendix D of Ref. [I29|. We stress the fact that the renormalization 
constant SZ W applies to both the W + and W~ field, i.e. the imaginary part of 5Z W is 
fixed by the renormalization condition and does not change sign when going from the W + 
to the W~ field. 

The renormalization conditions ()4.6|) have the solutions 



5/4 = E^(/4), 5nl = £f z (/4), 

2 2 

8Z Z A = — S T (0), SZ AZ = — j (/% 

A*z A^z 

5Z W = -Ef (/4), 5Z ZZ = -S^(/i 2 z ), 5Z AA = -E^», (4.7) 
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which require to calculate the self-energies for complex squared momenta. This would 
demand an analytic continuation of the 2-point functions entering the self-energies in the 
momentum variable to the unphysical Riemann sheet. In order to avoid this complication, 
we perform appropriate expansions about real arguments, 

S?(A4) = S t (K) + (Mw - M w )£f (M^) + 0(a 3 ), 
Ef*(/4) = Z ZZ (M 2 ) + 04 - M Z 2 )S^(M|) + 0(a% 

= -W(O) + ^f(Mi) - ^S^(O) + 0(a 2 ), (4.8) 

/i z [i z M z M z 

as done similarly in Appendix D of Ref. [I2H|- The 0(a 2 ) and 0(a 3 ) contributions result 
from products of terms Z w = 0(a), ^ zz = 0(a), ^ AZ = 0(a), (/4 - M&) = 0(a), 
and (/i| — M|) = 0(a) and are UV-fmite by construction at the one-loop level. 
By neglecting these 0(a 2 ) and 0(a 3 ) terms, we can replace ()4.7|) by 

5f4f = ^(M 2 ,) + 04 - M^)Sf (M 2 ■ 



^1 = S| Z (M Z 2 ) + (/4 - M Z 2 )S^(M Z 2 ), (4.9) 



and 



5Z ZA = (0), ^ = - A E ax (M 2 } + ( ^ - 1 ) *2 ZAj 

5Z W = -Sf (M^), = -S^(M Z 2 ). (4.10) 

We use these counterterms in our calculation. While the missing 0(a 2 ) terms in 5Z AZ 
do not influence our results, since the gauge-boson field renormalization constants drop 
out as there is no external gauge boson in the process under consideration, the missing 
(finite) 0(a 3 ) terms in the mass counterterms are beyond the accuracy of our calcula- 
tion. The counterterms ()4.9|) involve only functions that appear also in the usual on-shell 
renormalization scheme [EUEH], but consistently take into account the imaginary parts. 

When inserting the counterterms ()4.9|) and ()4.10|) into (J4.3)) , we can rewrite the renor- 
malized self-energies in the CMS as 

f™{k 2 ) = T^{k 2 ) - SM^ + (k 2 - M^)5Z W , 
t zz {k 2 ) = Y, zz {k 2 ) - 5M 2 + (k 2 - Mi)6Z zz , 
t AA {k 2 ) = Z AA {k 2 ) + k 2 5Z AA , 

t AZ (k 2 ) = Y, AZ (k 2 ) + k 2 Uz AZ + (k 2 - Ml)\bZ ZA (4.11) 



2 



with 



SM 2 , = £^(il4), SM 2 = £? Z (M 2 ), 



6Z ZA = ^S^(O), SZ AZ = -^L AZ {M 2 Z ) 



5Z W = -JT(M^), 5Z ZZ = -E'^(Mi), 5Z AA = -E'^(O). (4.12) 

Equations 1)4.11)1 with 1)4.12)1 have exactly the form of the renormalized self-energies in 
the usual on-shell scheme, but without taking the real part of the counterterms. While 
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in the on-shell scheme the self-energies are calculated in terms of the real renormalized 
masses M| and M^, in (|4.11j) and (|4.12|) the self-energies are to be calculated in terms 
of the complex masses /z| and /i^, although with real squared momenta. Note that this 
difference between usual on-shell and complex renormalization also changes the form of 
the IR divergence appearing in the W-field renormalization constant 5Zw- In the former 
scheme, it appears as logarithm lnm 7 of an infinitesimally small photon mass (or as the 
related 1/(4— D) pole in dimensional regularization); in the latter, the W width regularizes 
the singularity via lnT w . 

Owing to its definition (|4.1|) . the renormalization of the complex weak mixing angle is 
determined by 



5s w clSc^ cl (8n%{ 5$ 



(4.13) 



The electric charge is fixed in the on-shell scheme by requiring that there are no higher- 
order corrections to the ee'j vertex in the Thomson limit. In the CMS this condition reads 



e 2 c w /x z 



- = ;£ ,A4 (0) - f^^^. (4.14) 



Because of the presence of the complex masses in the loop integrals, the charge renormal- 
ization constant and thus the renormalized charge become complex. Since the imaginary 
part of the bare charge vanishes, the imaginary part of the charge renormalization con- 
stant is directly fixed by the imaginary part of self-energies. In a one-loop calculation, 
the imaginary part of the renormalized charge drops out in the corrections to the abso- 
lute square of the matrix element, because the charge factorizes from the lowest-order 
matrix element. Starting from the two-loop level, the imaginary part has to be taken into 
account. 

For a correct description of the resonances at the 0(a) level, we need the width 
including O(a) corrections. The width is defined via ()4.9|) . Using 5^ = M^ — /x^ and 
taking the imaginary part of (j4.9J) yields 

M w r w = Im{£^(M^)} - M w r w Re{Sf (M^)} + 0(a 3 ), (4.15) 

which can be iteratively solved for Tw- In 0(a 2 ), i.e. including first-order corrections to 
the width, the result is equivalent to the one obtained in the usual on-shell scheme. To this 
order the imaginary part of the self-energy is required in two-loop accuracy, but it can 
be more easily obtained by calculating the one-loop corrections to the decay processes 
W — > //' [ 29,. In our numerical calculation, we calculate the width from the decay 
processes including 0(a) corrections [H2|- 

In our calculation we consider external fermions only in the massless limit, in which 
these fermions are stable. Therefore, the corresponding on-shell self-energies do not in- 
volve any absorptive parts. Nevertheless they become complex via the complex renormal- 
ized weak mixing angle and the complex gauge-boson masses, and thus the field renor- 
malization constants 52f >cr of the fermion fields f a , defined by 

fS = (1 + ISZf^r, a = R,L, (4.16) 
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become complex. As in the case of 5Z\y, also these complex field renormalization constants 
apply both for fermions and antifermions, i.e. fields and antifields are not connected by 
complex conjugation anymore. Explicitly the SZf a are given by 

5Z La = -E f ' a (m 2 f ) - m) [E' f ' K (m 2 f ) + S // - L (mf) + 2£ ,/ ' s (mJ)] , (4.17) 

where a = R, L refers to the right- and left-handed components of the fermion self-energy 
£^(jo) following the conventions of Ref. [Hlj. Note that again no real part was taken in 
this relation in contrast to the usual on-shell renormalization. In contrast to SZw, there 
are soft IR divergences in the field renormalization constants 5Zf that are not regularized 
by finite widths but by the usual IR regulators such as m 7 . 

In the massless limit, the fermion-mass renormalization constant 5rrif tends to zero, 
and the quark-mixing matrix, if assumed to be different from the unit matrix, need not 
be renormalized. 

(ii) Complex renormalization for the Higgs sector and the top quark 

Since the top-quark and the Higgs boson do not appear in the lowest-order matrix 
elements for the processes under consideration, the corresponding mass and field coun- 
terterms are not needed in our calculation. Nevertheless, for completeness we define the 
corresponding renormalization constants here. 

For the Higgs boson the whole renormalization proceeds along the same lines as for 
the gauge bosons above. The complex Higgs mass squared 

/4 = M£ - iM H r H = M| - 5/4 (4. 18) 

is defined as the location of the zero in k 2 in the renormalized Higgs-boson self-energy 

± H (k 2 ) = Y> H (k 2 ) - 5/4 + (k 2 - fi 2 u )SZ H . (4.19) 

Fixing the Higgs field renormalization as above, we obtain the renormalization constants 



5& = X"(/4) = Z H (M 2 U ) + 04 - M*)E'*(m£) + 0(a% 

5Z H = -£'"(/4) = -X >H (M 2 ) + 0(a 2 ), (4.20) 

so that the renormalized Higgs-boson self-energy up to finite 0(a 2 ) terms can be written 

as 

t H {k 2 ) = Y, H {k 2 ) - 5M| + (k 2 - M£)6Z H . (4.21) 

with 

5Ml = £ H (M£), 5Z h = -E /H (M£). (4.22) 

For the top quark the renormalization procedure works analogously, i.e. we introduce 
the complex mass via 

H 2 t =m\ - im t r t , m tj0 = (H + Sfi t . (4.23) 
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The renormalized top-quark self-energy reads 

t\p) = V) + 5Ztp\ j)uj + + fs''V) + 5Z t , L ] fa. 



^(p 2 )-\{5Z t , K + 5Z t , h ) 



(4.24) 



in the conventions of Ref. [ I44j . Generalizing the on-shell renormalization conditions to 
complex renormalization as 



5Z t , a = -Vitf) 



cx = R,L, (4.25) 



fixes n\ as the location of the complex pole in p 2 in the top propagator. Expanding the 
self-energies appearing in the renormalization constants (|4.25|) about m\ and neglecting 
(UV-finite) terms of order 0(a 2 ), the renormalized top-quark self-energy can be expressed 
as 



£'(p) = £ t,R (p 2 ) + SZ t> n. i>u+ + E''V) + <^,l i>u 



+ fh 



5m t 



(4.26) 



with 



5m t 

sz t „ 



^ [E^Vt 2 ) + S*' L (m t 2 ) + 2S'- S (m t 2 ) 



S rt ' R (m 2 ) + £ / *< L (m t 2 ) + 2E 



a = R, L. 



(4.27) 



Finally, we complete the renormalization in the scalar sector. The field renormal- 
ization for the would-be Goldstone bosons can simply be set equal to the Higgs field 
renormalization constant 5Zh, which is sufficient to cancel all UV divergences in vertex 
functions. The tadpole counterterm 5t is again introduced to cancel explicitly occurring 
tadpole graphs, i.e. we set St = —t, where T H = it is the one-point vertex function for the 
Higgs boson at one loop. The gauge-fixing term need not be renormalized. 

In summary, in the CMS the usual renormalization conditions of the on-shell scheme 
can be used, but without taking any real parts. All parameters, in particular also the 
renormalization points, become complex. The calculation of self-energies with complex 
momentum arguments can be avoided by expanding these about real values. Performing 
these expansions appropriately and neglecting (UV-finite) higher-order terms, the renor- 
malized self-energies are obtained just as in the usual on-shell scheme with renormalization 
constants determined for real momenta. However, all the mass arguments are complex 
and no real parts must be taken. 

(in) Complex renormalization - background-field gauge 

Using the background-field method, the gauge-boson field renormalization constants 
can be determined in terms of the parameter renormalization in such a way that Ward 
identities possess the same form before and after renormalization [ 0U] . Real parameters 
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have to be substituted by the corresponding complex parameters everywhere when the 
complex renormalization is employed. The complex parameter renormalization is fixed 
as above in (jOjl . (jjX5j) , (l4"lij) . and (l4"2oT) . Note that S^ z (0) vanishes in the 

background-field gauge as a consequence of the background-field gauge invariance of the 
effective action, which in particular simplifies the charge renormalization constant ()4.14|) . 

Since the gauge-boson field renormalization constants drop out in the ^-matrix ele- 
ments without external gauge-boson fields, we can alternatively also use the definitions 
in (J4.1U)) in the calculation. 

(iv) Loop integrals with complex masses 

The consistent use of complex gauge-boson masses requires to use these also in the 
loop integrals. Thus, we need one-loop integrals with complex internal masses. The IR- 
singular integrals can be found in Ref. [HO]- Concerning the non-IR singular cases, we 
have analytically continued the results of Ref. [HB] for the 2-point and 3-point functions 11 , 
and the relevant results of Ref. [ [SUJ for the 4-point functions. We have checked all these 
results by independent direct calculation of the Feynman-parameter integrals. These 
results will be published elsewhere. 

5 Numerical results 

5.1 Input parameters and setup 

The numerical results are based on the same set of input parameters as in Ref. [137]: 



Gfj, 


= 1.16637 x KT 5 GeV- 2 , 


a(0) 


= 1/137.03599911, 


Oi s 


= 0.1187, 


M w 


= 80.425 GeV, 


M z 


= 91.1876 GeV, 


r z 


= 2.4952 GeV, 


M H 


= 115 GeV, 










m e 


= 0.51099892 MeV, 




= 105.658369 MeV, 


m T 


= 1.77699 GeV, 


m u 


= 66 MeV, 


m c 


= 1.2 GeV, 


m t 


= 178 GeV, 


m d 


= 66 MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV, 



which essentially follows Ref. [ 62j. For the top-quark mass mt we have taken the more 
recent value of Ref. [EH]- The masses of the light quarks are adjusted to reproduce 
the hadronic contribution to the photonic vacuum polarization of Ref. [El!- Since we 
parametrize the lowest-order cross section with the Fermi constant (G M scheme), i.e. 
we derive the electromagnetic coupling a according to ckg = V / 2G Ai M^(l — M^/M|)/7r, 
the results are practically independent of the masses of the light quarks. Moreover, 
this procedure absorbs the corrections proportional to m 2 /M^ in the fermion-W-boson 
couplings and the running of a{Q 2 ) from Q 2 = to the electroweak scale. In the relative 
radiative corrections, we use, however, a(0) as coupling parameter, which is the correct 
effective coupling for real photon emission. 

11 Note that the result of Ref. f° r the scalar two-point function is not valid in general for complex 
masses. In this case an extra rj function has to be added. The same comment applies to the results for 
the 2-point tensor integrals in Ref. [ I47| . 
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QCD corrections are treated in the "naive" approach of multiplying cross sections 
and partial decay rates by factors (1 + a s /^) per hadronically decaying W boson. The 
W-boson width T w is calculated from the above input including electroweak 0(a) and 
QCD corrections, yielding 

T w = 2.09269848 . . . GeV. (5.2) 

This procedure ensures that the effective branching ratios for the leptonic, semileptonic, 
and hadronic W decays, which result from the integration over the decay fermions, add 
up to 1. The value for the Z decay width T%, which is needed because of the Z resonance 
in the ISR convolution below the W-pair threshold, is taken from experiment [ :62 j • All 
other particles, including the top quark and the Higgs boson, are taken as stable. 

The setup differs from the one of Ref . [ [37] only in the event selection. In contrast to 
Ref. [E7J, where no phase-space cuts were applied at all, we now impose selection cuts 
in conjunction with a photon recombination procedure. In detail, we adopt the same 
procedure as in Ref. [l2*T|: 

1 . All bremsstrahlung photons within a cone of 5 degrees around the beams are treated 
as invisible, i.e. their momenta are disregarded when calculating angles, energies, 
and invariant masses. 

2. Next, the invariant masses M/ 7 of the photon with each of the charged final-state 
fermions are calculated. If the smallest M/ 7 is smaller than M rec = 25 GeV or if the 
energy of the photon is smaller than 1 GeV, the photon is combined with the charged 
final-state fermion that leads to the smallest Mf 7 , i.e. the momenta of the photon 
and the fermion are added and associated with the momentum of the fermion, and 
the photon is discarded 12 . 

3. Finally, all events are discarded in which one of the charged final-state fermions is 
within a cone of 10 degrees around the beams (after a possible recombination with 
a photon). No other cuts are applied. 

The presented results have been obtained with 10 8 events, using the subtraction 
method. All but the lowest-order predictions include naive QCD corrections and im- 
provements by ISR beyond 0(a), as described in Ref. [I21j. 

5.2 Results for differential cross sections 

In this section we consider results for various distributions at ^/s = 200 GeV and 
\/s = 500 GeV in the setup described above. For reference, in Table [D we provide the 
corresponding integrated cross sections for the final states v t t + fi^v^, udfi^u^, and udsc 
in various approximations for different CM energies ^Js. The numbers in parentheses 
represent the uncertainties from Monte Carlo integration in the last digits of the predic- 
tions. Columns two and three in each table contain the two versions of the lowest-order 
cross section for the full e + e~ — > 4/ processes corresponding to the different treatments of 
finite-width effects as provided by the "fixed-width scheme" (FW) and the complex-mass 
scheme (CMS). In the FW scheme the finite constant width, and thus the complex mass, 

12 Except for the 1 GeV cut, the described cut and recombination procedure coincides with the one used 
in the first two papers of Ref. [ I2()j . 
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v/i/GeV Born(FW) 


Born(CMS) 


DPA 


ee4f 


e + e — > v t t + n 


200 211.52(3) 


211.40(3) 


191.98(3) 


191.18(3) 




[-0.06%] 


[-9.24(1)%] 


[-9.57(1)%] 


500 62.17(1) 


62.14(1) 


65.48(2) 


64.96(2) 




[-0.05%] 


[+5.32(1)%] 


[+4.54(1)%] 


e + e — > ud/i Vfj, 


200 628.72(9) 


628.37(9) 


591.55(9) 


589.11(10) 




[-0.06%] 


[-5.91(1)%] 


[-6.25(1)%] 


500 180.83(4) 


180.73(4) 


197.87(5) 


196.20(5) 




[-0.06%] 


[+9.42(1)%] 


[+8.56(2)%] 


e + e — > udsc 


200 1868.9(3) 


1867.8(3) 


1822.6(3) 


1815.3(3) 




[-0.06%] 


[-2.48(1)%] 


[-2.81(1)%] 


500 526.1(1) 


525.8(1) 


597.1(2) 


591.9(2) 




[-0.06%] 


[+13.50(2)%] 


[+12.58(2)%] 



Table 1: Integrated cross sections in fb for e + e~ — ► u t t + ijTv^ ud/i"^, and udsc in Born 
approximation (in the fixed-width and complex- mass schemes), in DPA, and using the 
full 0(a) correction (ee4f); all but the Born cross sections include higher-order ISR and 
(if relevant) naive QCD corrections. 

is only inserted into the propagators. The relative difference cr Born (CMS)/0"Born(FW) — 1 
of the schemes in lowest order is given by the numbers in square brackets in the third 
columns. We have not given an error on this difference, because the two Born predictions 
are strongly correlated. The last-but-one columns show the DPA of RacoonWW which 
also includes some effects beyond DPA, as described in Ref. [|2T]; the numbers in square 
brackets are defined as £dpa = o"DPA/ cr Bom(FW) — 1. We normalize odpa to crBorn(FW), 
because the lowest-order part of the DPA is per default evaluated in the FW scheme in 
RacoonWW. Finally, the last columns (ee4f) contain the full one-loop corrections to 
e + e~ — > 4/; the numbers in square brackets are defined as 5 ee 4f = <7 ee 4f/o'Born(CMS) _ 1- 
Here we normalize to 0"B O rn(CMS), because the full e + e~ — > 4/ calculation is consistently 
performed in the CMS. For y/s = 500 GeV the difference between the predictions based 
on the full 0(a) corrections and on the DPA slightly increases by ~ 0.1% with respect to 
the results without cuts [EZ|- This tendency can be attributed to the fact that at high 
energies the cross section is dominated by forward-scattered nearly on-shell W-bosons. 
Such events are discarded by the cuts thus reducing the contribution of on-shell W-boson 
pairs and worsening the quality of the DPA. 

For the differential distributions, we focus on the semileptonic process e + e~ — > udfi'u^ 
in the following. The respective results for the final states u t t + ^u^ and udsc look similar, 
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up to an offset resulting from the QCD corrections; in particular, the difference between 
the DPA and the full 0{a) calculation are almost identical. We always display the 
lowest-order prediction (Born) and the result of the full one- loop calculation (ee4f) in the 
upper row of each figure. The relative corrections (in per cent) in the DPA approach 
(DPA), <5dpa = do"DPA/do"Bom(FW) — 1, and in the full one-loop calculation (ee4f), <5 ee4 f = 
dcr ee 4f/dcrBom(CMS) — 1, are shown in the lower rows of the figures. These additionally 
include an inset depicting the relative difference between the full one-loop and the DPA 
calculation with respect to the DPA calculation, A = dcx ee 4f /d<7DPA — 1- We define all 
angles in the laboratory system, which is the CM system of the initial state. The momenta 
of the W + and W~ bosons are defined as 



k + = ki + k 2 , k- = k 3 + k 4 , (5.3) 

respectively, after a possible photon recombination. From these momenta the invariant 
masses of the virtual W bosons and their angles are calculated. 

The invariant-mass distributions for the W + and W~ bosons are shown in Figures El 
and From the plots for the relative corrections, it can be seen that the full one-loop 
corrections are smaller than the DPA corrections for invariant masses bigger than M\y 
and vice versa for invariant masses smaller than M^. If neglected, this effect will give rise 
to a small shift in the direct reconstruction of the W-boson mass. The distribution in the 
cosine of the W + production angle 6\v+ is shown in Figure ED While at LEP energies there 
is hardly any distortion of the shape induced by corrections beyond DPA, at 500 GeV the 
difference of the corrections in DPA and the complete 0(a) corrections rises from —1% to 
about —2.5% with increasing scattering angle. Note that such a distortion of the shape of 
the angular distribution can be a signal for anomalous triple gauge-boson couplings. The 
angular dependence of A is even more pronounced in the distribution in the decay angle 
&w-fi- presented in Figure [T3J Note, however, that as a general feature the cross section 
is smallest where the corrections beyond DPA are largest. The distribution in the energy 
E^- of the muon can be found in Figure Again, for -^i = 200 GeV the corrections 
beyond DPA hardly depend on the muon energy in the interval 20 GeV < < 80 GeV, 
where two resonant W bosons are kinematically possible. Outside this interval, the DPA 
runs out of control, and the difference to the full calculation becomes big; but there the 
cross section is very small. 

Finally, we consider the distributions in azimuthal angles that were also discussed in 
Ref . [ |^j including corrections in DPA. In Figure E3 we show the distributions in the 
azimuthal decay angle 0w+ °f t ne W + boson, i.e. the angle between the decay plane of 
the W + and the plane of W-pair production, 

(k+ x p+)(k+ x ki) 

COS0 W + - 



|k + x p+||k + x ki| 
sgn(sin0 w +) = sgn{k + ■ [(k + x p + ) x (k+ x ki)]} . (5.4) 

Here, in particular at 200 GeV the difference between DPA and the full 0(a) corrections 
is approximately proportional to sin0w+ phis some constant offset. In the DPA, it can be 
deduced from Appendix A of Ref. [ HH] , that the contributions of the imaginary parts of the 
one-loop coefficient functions always involve a factor sin 0w+ o r a factor sin 0w- together 
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Figure 10: Distribution in the invariant mass of the W + boson (upper row) and the 
corresponding corrections (lower row) at yfs = 200 GeV (l.h.s.) and y/s = 500 GeV 
(r.h.s.) for e + e~ — > ud/i^z/^. The inset plot shows the difference between the full 0(a) 
corrections and those in DPA. 
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d^T [(Osv\ = 200 GeV cofcr [mv] = 500 GeV 




75 76 77 78 79 80 81 82 83 84 85 75 76 77 78 79 80 81 82 83 84 85 

M w - [GeV] M w - [GeV] 

Figure 11: Distribution in the invariant mass of the W~ boson (upper row) and the 
corresponding corrections (lower row) at yfs = 200 GeV (l.h.s.) and y/s = 500 GeV 
(r.h.s.) for e + e~ — > ud/i^z/^. The inset plot shows the difference between the full 0(a) 
corrections and those in DPA. 
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dcose w + [fb] v / s = 200GeV 



dcose w+ N v^ = 500GeV 



1000 - 




1000 - 




cos # w + 



COS # w + 



Figure 12: Distribution in the cosine of the W + production angle with respect to the 
e + beam (upper row) and the corresponding corrections (lower row) at ^/s = 200 GeV 
(l.h.s.) and y/s = 500 GeV (r.h.s.) for e + e~ — > vAji~u^. The inset plot shows the 
difference between the full 0(a) corrections and those in DPA. 
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dcosC-„- ^ v^ = 200GeV O^:^ v^ = 500GeV 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

Figure 13: Distribution in the cosine of the \i~ decay angle with respect to the W~ 
direction (upper row) and the corresponding corrections (lower row) at \fs = 200 GeV 
(l.h.s.) and y/s = 500 GeV (r.h.s.) for e + e~ — > vAji~u^. The inset plot shows the 
difference between the full 0(a) corrections and those in DPA. 
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do- f fb 1 dcr I" fb 1 

dE~ [Gev\ yfe = 200 GeV dE~ [Gev\ ^ = 500 GeV 




10 20 30 40 50 60 70 80 90 100 50 100 150 200 250 



Ep- [GeV] E tl - [GeV] 

Figure 14: Distribution in the energy of the pT (upper row) and the corresponding 
corrections (lower row) at ^fs = 200 GeV (l.h.s.) and y/s = 500 GeV (r.h.s.) for 
e + e~ — > ud/i^z/^. The inset plot shows the difference between the full 0(a) corrections 
and those in DPA 
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with symmetric functions in these angles. Contributions of real parts, on the other hand, 
are symmetric in these angles. In order to illustrate this feature, we have included an 
extra curve in the plots for the relative corrections to this distribution labelled "DPA 
(real)". In this curve, we have switched off all imaginary parts in the DPA calculation, 
apart from the finite width in the resonant propagators. As expected, these results are 
symmetric in the angle 0w+ about 180°. The contribution proportional to sin0 w + in 
the difference of the full 0(a) calculation with respect to the DPA (real) is even larger 
than the corresponding difference with respect to the complete DPA. These properties 
of the DPA suggest, that also the sinusoidal dependence on W + m the latter difference 
results from imaginary parts. A possible source could be the missing imaginary parts of 
the counterterms in DPA. 

The distributions in the angle between the two planes spanned by the momenta 
of the two fermion pairs in which the W bosons decay, i.e. (note that k + = — k for 
non-photonic events) 

(k+ x ki)(-k_ x k 3 ) 

COS( P = Ti 1 ii i TT> 

|k + x ki||— k_ x k 3 | 

sgn(sin^) = sgn {k + ■ [(k + x ki) x (— k_ x k 3 )]} , (5.5) 

are presented in Figure [TBI The large corrections for angles <p near 0° or 180°, i.e. if the 
two decay planes coincide, result from the suppression of hard photonic corrections [ |^] . 
The corrections beyond DPA do hardly depend on 0. 

6 Conclusions 

We have presented technical and conceptual details as well as further numerical results 
of a calculation of the complete electroweak 0(a) corrections to the charged-current four- 
fermion production processes e + e~ — ► v T r + "fJTVp, ud/i'z/^, and udsc. 

In particular, we have described methods how the O(10 3 ) occurring different spinor 
structures can be algebraically reduced to a few simple standard structures. The pre- 
sented algorithms, which shorten the analytical results considerably and thereby render 
the resulting computer code relatively short, should be useful for other future calculations. 

Moreover, a concept for consistently performing one-loop calculations with complex 
masses for unstable particles is presented. Technically the complex masses are introduced 
via an appropriate complex renormalization prescription at the level of the Lagrangian, 
so that the usual machinery of perturbation theory (Feynman rules etc.) can be used 
to organize the calculation. Since the theory is not changed at all, there is no danger of 
double counting terms. The complex mass parameters, which also enter the weak mixing 
angle and coupling constants, can be viewed as an analytical continuation of the real 
masses. Consequently, all algebraic relations that follow from gauge invariance (Slavnov- 
Taylor identities, Ward identities, cancellation of gauge-parameter dependences) are valid 
in spite of the complex masses which incorporate finite-width effects. The price to pay for 
this achievement is that all one-loop integrals have to be performed with complex mass 
parameters. 

Finally, we have completed our discussion of numerical results, which was started in 
Ref. [EZ1 by considering total cross sections, by showing the effects of the complete 0(a) 
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Figure 15: Distribution in the azimuthal decay angle of the W + (upper row) and the 
corresponding corrections (lower row) at = 200 GeV (l.h.s.) and -/s = 500 GeV (r.h.s.) 
for e + e~ — > ud/i~z/ M . The inset plot shows the difference of the full 0(a) corrections to 
those in DPA and the larger difference to those in DPA (real), which is calculated without 
imaginary parts. 
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Figure 16: Distribution in the azimuthal angle (upper row) and the corresponding 
corrections (lower row) at -^i = 200 GeV (l.h.s.) and t/s = 500 GeV (r.h.s.) for e + e~ — > 
ud/^~z/ M . The inset plot shows the difference between the full 0(a) corrections and those 
in DPA. 
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corrections to various differential cross sections of physical interest. In particular, we have 
considered differential cross sections sensitive to imaginary parts. Our results have also 
been compared to predictions based on the double-pole approximation, revealing that 
the latter approximation is not sufficient to fully exploit the potential of a future linear 
collider in an analysis of W-boson pairs at high energies. Specifically, at (and above) a 
CM energy of 500 GeV the corrections beyond DPA induce a non-negligible distortion in 
angular distributions that could be misinterpreted as signal for anomalous triple gauge- 
boson couplings if not taken into account. The remaining theoretical uncertainties for 
the total cross section have been discussed in Ref. [EHj. For differential distributions the 
uncertainties from QCD effects are even more relevant. 

In this paper we have presented methods that were successfully used for the first com- 
plete 0(a) calculation of electroweak corrections to a process with six external particles 
and involving unstable particles in intermediate states. Processes of this type will become 
more and more important in the future. Our methods are not specifically adapted to the 
considered process and therefore should be helpful in precision calculations for similar 
processes at future colliders. 
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Appendix 

Interpretation of the reduction algorithm of Section 13.21 in the Weyl— van der 
Waerden spinor formalism 

In this appendix we inspect the reduction algorithm of Section l3~2l within the Weyl- 
van der Waerden (WvdW) spinor formalism, which is particularly simple for massless 
particles. We consistently employ the conventions of Ref. [EH] for the spinor method and 
highlight only those features that are crucial for our manipulations. 

In this formalism, Dirac spinors and momenta of massless fermions are described by 
Weyl spinors pa, Pa = (Pa)* which are related to the four-momentum p^ according to 



where a^ B denote the unit matrix (/x = 0) and the Pauli matrices (/x = 1,2,3). The 
matrices cr^ 1 also appear as the non- vanishing blocks in the Dirac matrices 7^. Thus, 
Dirac spinor chains translate into chains involving Weyl spinors and cr^ 1 matrices. Coupled 
chains can be separated by using the relation 



P» a AB = PaPb, 



(A.l) 



a AB a v,CD - 2e AC 6 BD, 
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where €ab is the totally antisymmetric tensor in two dimensions and e A ^ is its complex 
conjugate. In this way, all Weyl spinors get contracted in so-called spinor products 

(PiPj) = £ AB Pi,APj,B (A. 3) 

or their complex conjugates. Using (|A.2jl . Lorentz products can also be translated into 
spinor products. We illustrate the above considerations by the simplest examples: 

(PiPj) = \(PiPj)(PiPj)\ 



fa .. = (pkPi){PkPjY 



V 



fa .. = (PkPi)*(PkPj)- (A.4) 



ij 



Dirac chains of the form ()3.3|) involving no additional e^ upa can thus be directly reduced 
to a product of (at least four) spinor products. The totally antisymmetric tensor e^ upa is 
translated into WvdW objects according to 



z pVpK 



W H u,BC k p,DA _ fi p,BC « n v,DA\ ( a k\ 

A \ a AB a a CD G G AB a a CD a ) > l A - D J 



so that contractions with e pupa can also be expressed in terms of WvdW spinor products. 
The two-dimensionality of the WvdW spinor space is rather restrictive and leads to the 
identity 

(PiPj) (PkPi) + (PiPk) (PiPj) + (PiPi) (PjPk) = 0, (A.6) 

which frequently admits simplifications in complicated expressions. In the cases relevant 
for the considered processes, these identities allow to combine the sums that result from 
the application of (jA.5|) . so that also all spinorial expressions involving contractions with 
e pu P a can k e rec i ucec i to a product of spinor products. 

The identities (jA.5|) and (jA.6|) also allow to write the combinations A^^i °f f)3.21|) as 
simple factors of spinor products: 

A tjkl~ = 2 (P^j ) &kPl ) (PiPk ) * (PjPl ) * , 

A tjkl = \{PtPl) l <PkPj){PiPjY{PkPl)\ 
1 

A ijkt = ^(PiPk)(PjPl)(PiPl)*(PjPk)*, 

Kli = {At-itT- (A.7) 

When expressing the Af^f in terms of spinor products, the relations (|3.3U|) - (|3.33|) 
become trivial. The relations (|A.7|) and the fact that in the WvdW formalism all spinor 
chains can be expressed in terms of products of spinor products explains why in the 
formalism of Section 13.21 all spinor chains could be reduced to a simple product of Afj^f 
and scalar products. The results of ()3.34|) read in the WvdW formalism 



12 



^plu^p 



34 = S2(p 2 P4)(pm)*S a+ 6 T+ + 32(p 2 P4)*(PlP3)5a-5r- 
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+ 1 



hi Vl 12 \ivYi\ , JtmTptJ kk = 32(^5) (p 4 p 6 ) 
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5(i 



X (<J<T+(P2P3><PlPs)* + <k-(PlP3)(P2P3)*) , 
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foYYYl" 12 7^7p7« „„ 7/, Kfi = 64(pip 3 )(p 3 p 5 ) (p 2 Pi)* (PsPe Y 



r >(i 



-ie^> liCT [ 7/ ,J 12 [7,J 34 [7pJ 56 = -2(piP2)(j>3P 5 ){pm)*{pip e )* • (A.8) 
Finally, the relevant SMEs ()3.26|) read in the WvdW formalism 

= {pm){p2P-iY{pm){piPiY{p\Pb){pm)\ 
M + ~~ = (pmY (P2P3) ipm) (pmY (pips) (pmY ■ ( A - 9 ) 

Of course, it would have been possible to translate each spinor chain into WvdW spinor 
products directly after performing the loop integration, i.e. skipping the steps described 
in Section l3~21 However, it would have been hard to perform all the manipulations with 
WvdW objects made there in terms of Lorentz products. For instance, it is practically 
impossible to make full use of momentum conservation in very involved expressions in 
terms of WvdW spinor products. It turned out in many examples given in the literature 
that no algorithmic recipe has been found yet yielding the most compact expressions 
in complicated amplitudes. Moreover, the generalization to massive fermions is more 
complicated in the WvdW formalism than in the method described in Section 13.21 In fact 
a variant of this method has already been used in Ref. [03] for the process e + e~ — > ttH. 
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